The scalar bi-spectrum during preheating in single field inflationary models 
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In single field inflationary models, preheating refers to the phase that immediately follows in- 
flation, but precedes the epoch of reheating. During this phase, the inflaton typically oscillates 
at the bottom of its potential and gradually transfers its energy to radiation. At the same time, 
the amplitude of the fields coupled to the inflaton may undergo parametric resonance and, as a 
consequence, explosive particle production can take place. A priori, these phenomena could lead 
to an amplification of the super-Hubble scale curvature perturbations which, in turn, would modify 
the standard inflationary predictions. However, remarkably, it has been shown that, although the 
Mukhanov-Sasaki variable does undergo narrow parametric instability during preheating, the am- 
plitude of the corresponding super-Hubble curvature perturbations remain constant. Therefore, in 
single field models, metric preheating does not affect the power spectrum of the large scale pertur- 
bations. In this article, we investigate the corresponding effect on the scalar bi-spectrum. Using 
the Maldacena's formalism, we analytically show that, for modes of cosmological interest, the con- 
tributions to the scalar bi-spectrum as the curvature perturbations evolve on super-Hubble scales 
during preheating is completely negligible. Specifically, we illustrate that, certain terms in the third 
order action governing the curvature perturbations which may naively be expected to contribute 
significantly are exactly canceled by other contributions to the bi-spectrum. We corroborate selected 
analytical results by numerical investigations. We conclude with a brief discussion of the results we 
have obtained. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Inflation is currently considered the most attractive 
paradigm for explaining the extent of homogeneity of 
the observable universe. In addition, the paradigm also 
provides a well motivated causal mechanism for the gen- 
eration of perturbations in the early universe (see any 
of the following texts [1-7] or one of the following re- 
views [8-20]). In the standard picture, at the end of the 
inflationary epoch, the inflaton, which is coupled to the 
other fields of the standard model, decays into relativistic 
particles thereby transferring its energy to radiation (see, 
for instance, Refs. [12, 18] and also Refs. [21-24]). These 
decay products are then expected to thermalize [25] in 
order for the radiation dominated epoch corresponding 
to the conventional hot big bang model to start. 

In many models, inflation is terminated when the 
scalar field has rolled down close to a minimum of the 
potential. Thereafter, the scalar field usually oscillates 
at the bottom of the potential with an ever decreasing 
amplitude. There exists a period during this regime, im- 
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mediately after inflation but prior to the epoch of reheat- 
ing, a phase that is often referred to as preheating [26- 
28]. During this brief phase, as in the inflationary era, 
the scalar field continues to remain the dominant source 
that drives the expansion of the universe. Though the 
modes of cosmological interest (corresponding to comov- 
ing wavenumbers k such that, say, 10 -4 < k < 1 Mpc -1 ) 
are well outside the Hubble radius during this phase, 
the conventional super-Hubble solutions to the curvature 
perturbations that are applicable during inflation do not 
a priori hold at this stage. In fact, careful analysis is 
required to evolve these modes during the phase of pre- 
heating. However, despite the subtle effects that need 
to be accounted for, it can be shown that, in single field 
inflationary models, the amplitude of the curvature per- 
turbations and, hence, the amplitude as well as the shape 
of the scalar power spectrum associated with the scales of 
cosmological interest remain unaffected by the process of 
preheating (for the original effort, see Ref. [29]; for more 
recent discussions, see Refs. [30-32]). 

Over the last decade, there has been a tremendous the- 
oretical interest in understanding the extent of the non- 
Gaussianities that are generated during inflation [33-57]. 
Simultaneously, there has been a constant effort to ar- 
rive at increasingly tighter constraints on the dimension- 
less non-Gaussianity parameter / NL that is often used 
to characterize the amplitude of the reduced scalar bi- 
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spectrum (viz. a suitable ratio of the scalar bi-spectrum 
to the corresponding power spectrum) from the available 
Cosmic Microwave Background (CMB) data [58-76]. For 
instance, it has been theoretically established that slow 
roll inflation driven by the canonical scalar field typically 
leads to rather small values of / NL (of the order of the 
first slow roll parameter) [33-35, 37]. In contrast, though 
a Gaussian primordial perturbation lies well within 2-cr, 
the mean values of / NL from the CMB observations seem 
to indicate a significant amount of non-Gaussianity [58- 
76] . Such large levels for the parameter / NL can be gener- 
ated when one considers non-canonical scalar fields [38- 
42] or when there exist deviations from slow roll infla- 
tion [44-52]. 

We mentioned above that, in single field models, the 
epoch of preheating does not affect the curvature pertur- 
bations and the scalar power spectrum generated during 
inflation on cosmologically relevant scales. Note that, 
the scalar power spectrum is essentially determined by 
the amplitude of the curvature perturbation. Whereas, 
as we shall discuss, the scalar bi-spectrum generated dur- 
ing inflation involves integrals over the curvature pertur- 
bations as well as the slow roll parameters [37-42, 44- 
52]. If indeed deviations from slow roll inflation can re- 
sult in high levels of non-Gaussianity, then, naively, one 
may imagine that the termination of inflation and the 
regime of preheating — both of which involve large val- 
ues for the slow roll parameters — can also lead to large 
non-Gaussianitics. In other words, it may seem that the 
epoch of preheating can contribute significantly to the 
scalar bi-spectrum. In this work, we shall investigate the 
contributions to the scalar bi-spectrum during preheating 
in single field inflationary models. Remarkably, though 
the epoch of preheating actually amplifies specific contri- 
butions to the bi-spectrum, as we shall illustrate, certain 
cancellations arise that leave the total bi-spectrum gen- 
erated during inflation virtually unaltered. 

This paper is organized as follows. In the following 
section, we shall highlight the essential aspects of pre- 
heating in single field inflationary models. In particular, 
we shall discuss the behavior of the scalar field as well as 
the scale factor, when a canonical scalar field is oscillating 
at the bottom of an inflationary potential which behaves 
quadratically near its minimum. We shall also discuss a 
few important points concerning the evolution of the cur- 
vature perturbation on supcr-Hubble scales during pre- 
heating. In Sec. Ill, after defining the bi-spectrum, we 
shall outline the various contributions to the bi-spectrum 
in the Maldacena formalism. In Sec. IV, we shall evalu- 
ate the contributions to the bi-spectrum for super- Hubble 
modes as the scalar field oscillates in the quadratic po- 
tential, and show that the total contribution during this 
epoch proves to be insignificant for these modes. We 
shall also support certain analytical results with the cor- 
responding numerical computations. Finally, in Sec. V, 
we shall conclude with a brief summary and outlook. 

A couple of words on our notation are in order at this 
stage of our discussion. We shall work with units such 



that c = tl = 1 and shall set = (8 7TG)" 1 . An over- 
dot and an overprime shall denote differentiation with 
respect to the cosmic time t and the conformal time 77, 
respectively. Also, N shall represent the number of e- 
folds. Moreover, double angular brackets shall denote 
averaging over the oscillations during preheating. 



II. BEHAVIOR OF THE BACKGROUND AND 
THE LARGE SCALE PERTURBATIONS 
DURING PREHEATING 

In this section, we shall discuss the behavior of the 
background and the evolution of the curvature perturba- 
tion on super-Hubble scales during preheating. We shall 
consider a model involving the canonical scalar field and 
assume that the inflationary potential behaves quadrati- 
cally around its minimum. 



A. Background evolution about a quadratic 
minimum 

Consider a canonical scalar field </> that is governed by 
the quadratic potential V(<p) = m 2 cf> 2 /2 near its min- 
imum. It is well known that, in such cases, slow roll 
inflation can be realized if the field starts sufficiently far 
away from the minimum, with suitably small values for 
its velocity [1-20]. Provided the initial conditions fall 
in the basin of the inflationary attractor, the number of 
c-folds of inflation achieved largely depends only on the 
initial value of the field, and inflation ends as the field 
nears the bottom of the potential. In fact, according 
to the slow roll approximation, in an inflationary po- 
tential that consists of no terms other than the above- 
mentioned quadratic one, inflation gets terminated as 
the field crosses O = v2 M pl . Thereafter, the scalar 
field oscillates about the minimum with a constantly de- 
creasing amplitude because of the friction caused due to 
the expansion. These behavior are clearly evident from 
Fig. 1, where we have plotted the evolution of the scalar 
field and the first slow roll parameter e 1; arrived at nu- 
merically, both during and immediately after inflation for 
the quadratic potential. We should emphasize here that 
focusing on single field models can be considered to be 
essentially equivalent to assuming that the coupling of 
the inflaton to other fields is suitably weak during pre- 
heating. The weak coupling will allow the condensate to 
live sufficiently long for a few oscillations to take place 
about the minimum of the inflaton potential. 

Recall that, the first slow roll parameter is given by 



H 2 2H 2 M 2 ' w 
pi 

where H = a/a is the Hubble parameter, with a(t) being 
the scale factor associated with the Friedmann-Lemaitre- 
Robertson- Walker line-element. The second slow roll pa- 
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FIG. 1: The behavior of the scalar field (left panel) and the evolution of the first slow roll parameter ei (right panel) during 
the epochs of inflation and preheating have been plotted as a function of the number of e-folds for the case of the archetypical 
chaotic inflationary model described by the quadratic potential. The blue curves denote the numerical results, while the dotted 
red curves in the insets represent the analytical results given by Eqs. (4) and (6) that are applicable during preheating. The 
analytical results evidently match the numerical ones quite well. Note that, for the choice parameters and initial conditions 
that we have worked with, ei turns unity at the e-fold of N c ~ 28.3, indicating the termination of inflation at the point. The 
fact that the field oscillates with a smaller and smaller amplitude once inflation has ended is clear from the inset (in the figure 
on the left panel). We should mention that we have worked with a smaller range of e-folds just for convenience. 



rameter €2 is denned in terms of the first one as follows: scalar field satisfies the differential equation 
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Since the field is oscillating at the bottom of the poten- 
tial post inflation, <j) = at the 'turning points' and, 
hence, t\ = at such instances. Also, when the field is 
at the bottom of the potential, V((f> = 0) = so that 
QH 2 = 4> 2 , corresponding to ei = 3. Hence, we can 
expect ei to oscillate between these two extreme values. 
Moreover, the above expression for the second slow roll 
parameter in terms of the first suggests that £2 will van- 
ish whenever t\ reaches the maximum value (i.e. at the 
bottom of the potential wherein l\ = 0), and that it 
will diverge at the 'turning points' wherein e\ vanishes. 
These behavior too are indeed reflected in the plot of ei 
in Fig. 1 and, in Fig. 2, where we have plotted 62- 

Let us now try to arrive at the complete behavior of the 
background scalar field analytically. During the phase 
of preheating, one finds that the period of the oscilla- 
tions (characterized by the inverse mass in the case of the 
quadratic potential of our interest) is much smaller than 
the time scales associated with the expansion, i.e. the in- 
verse of the Hubble parameter H . In such a situation, 
to understand the effects of the scalar field on the scale 
factor, one can average over the oscillations and make 
use of the averaged energy density of the scalar field to 
solve the first Friedmann equation. One finds that, in 
the quadratic potential of our interest, the expansion be- 
haves as in a matter dominated epoch, with the scale 
factor growing as a(t) cx i 2 / 3 , so that the Hubble param- 
eter behaves as H = 2/(3 1) [12, 18, 21, 22]. 

Therefore, at the time of preheating, the oscillating 
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The solution to this differential equation can be immedi- 
ately written down to be 



-7^7— = — sin (mt + A), 
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where a is a dimensionless constant that we shall soon 
determine, while A is an arbitrary phase chosen suitably 
to match the transition from inflation to the matter dom- 
inated era. The 'velocity' of the field is then given by 
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cos (m t + A) sin (m t + A) 

mt 



— cos (mt + A), 



(5) 



where, in arriving at the second expression, for the sake 
of consistency (i.e. in having made use of the averaged 
energy density in the first Friedmann equation to arrive 
at the scale factor), we have ignored the second term in- 
volving t~ 2 . Upon using the above expressions for r/> and 
(f) and the fact that H = 2/(3 1) in the first Friedmann 
equation, we obtain that a = ^8/3. Under these condi- 
tions, we find that the first slow roll parameter simplifies 
to 



ei ~ 3 cos 2 (mt + A) 



(0) 



which, upon averaging, reduces to 3/2, as is expected in 
a matter dominated epoch. It is represented in Fig. 1 
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B. Evolution of the perturbations 

As is well known, the scalar perturbations are essen- 
tially described by the curvature perturbation, say, 1Z. 
The scalar power spectrum V s (k) can be defined through 
the two point correlation function of the Fourier modes 
of the curvature perturbation as follows [8-11, 16-20] 



(2^ 
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V s (k) S {3) (fe+p). 



(8) 



The quantum operator TZk in turn is decomposed based 
on the modes fk which are governed by the differential 
equation 



FIG. 2: The behavior of the second slow roll parameter ti 
immediately after the termination of inflation has been plot- 
ted as a function of e-folds. As in the previous figure, the 
blue curve represents the numerical result, while the dashed 
red curve denotes the analytical result during preheating 
[viz. Eq. (7)]. Upon comparing this plot with the earlier plot 
of £i , it is clear that £2 diverges exactly at the turning points 
where ei vanishes, while £2 itself vanishes whenever the field 
is at the bottom of the potential at which point £i attains its 
maximum value. 



(right panel). On using the above result for e\ in the 
definition (2), the second slow roll parameter can be ob- 
tained to be 



fk+2-fL + k 2 h=0, 



(9) 



e 2 {t) ~ — 3mt tan (mt + A) 



(7) 



which is illustrated in Fig. 2. 

During preheating, we can write a(t) = a c (t/t c ) 2 / 3 , 
where t e and a e denote the cosmic time and the scale 
factor at the end of inflation. We should mention here 
that, in addition to the phase A, one requires the value 
of t c in order to match the above analytical results for 
the scalar field and the slow roll parameters with the 
numerical results. After setting t c — 7 [2/(3 H c )], where 
H c is the value of the Hubble parameter at the end of 
inflation, we have chosen the quantity 7 and the phase A 
suitably so as to match the analytical expressions with 
the numerical results. It is clear from Figs. 1 and 2 that 
the agreement between the numerical and the analytical 
results is indeed very good. We should mention here 
that, for the results to match, we seem to require a 7 
that is slightly larger than unity. Actually, for the values 
of the parameters that we have worked with, we find that 
we need to choose 7 to be about 1.18 for the analytical 
results to match the numerical ones. The fact that 7 is 
not strictly unity need not come as a surprise. After all, 
some time is bound to elapse as the universe makes the 
transition from an inflationary epoch to the behavior as 
in a matter dominated era. 



where z = y/2 1\ M pl a. The scalar power spectrum can 
be expressed in terms of the modes fk as 
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(10) 



with the assumption that the spectrum is evaluated at 
suitably late times when the modes are sufficiently out- 
side the Hubble radius during the era of inflation [1- 
4, 6, 8-11, 16-20]. 

It often proves to be convenient to work in terms of the 
so-called Mukhanov-Sasaki variable Vk, which is defined 
as = z fk- In terms of the variable Vk , the above 
equation of motion for fk reduces to the following simple 
form: 
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v k 



k z - — 

z 



v k = 0. 



The initial conditions on the perturbations are imposed 
when the modes are well inside the Hubble radius dur- 
ing inflation. The modes are usually chosen to be in 
the Bunch-Davies vacuum, which amounts to demanding 
that the Mukhanov-Sasaki variable Vk reduces to follow- 
ing Minkowski-likc positive frequency mode in the sub- 
Hubblc limit: 



lim Vk 

k I (a H)—}oo 



1 



-i krj 



(12) 



As we shall see, the scalar bi-spectrum shall involve in- 
tegrals over the modes fk and its derivative f' k as well as 
the slow roll parameters ei, £2 and the derivative e 2 . So, 
in order to analyze the effects on the bi-spectrum due to 
preheating, it becomes imperative that, in addition to the 
behavior of the slow roll parameters, we also understand 
the evolution of the mode fk and its derivative during 
this epoch. We have already studied the behavior of the 
first two slow roll parameters in the previous sub-section. 
Therefore, our immediate aim will be to understand the 
evolution of the curvature perturbations for scales of cos- 
mological interest during the preheating phase. 
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Since the modes of cosmological interest are well out- 
side the Hubble radius [i.e. k/(aH) <C 1] at late times, we 
need to arrive at the super-Hubble solution for the mode 
fk or, equivalently, the Mukhanov-Sasaki variable Vk- In 
a slow roll inflationary regime, i.e. when (£1,62,63) <C I, 
the effective potential z" / z that governs the evolution of 
Vk [cf. Eq. (11)] can be written as 



a 2 H 2 [2 + 0{e 1 ,6 2 ,e 3 )}~2a 2 H 2 . 



(13) 



Due to this reason, during slow roll inflation, the supcr- 
Hubble condition k/(aH) <C 1 amounts to neglecting 
the k 2 term with respect to the effective potential z" Jz 
in the differential equation (11). In such a case, it is 
straightforward to show that the super-Hubble solution 
to Vk up to the order k 2 can be expressed as follows [1- 
4, 6, 8-11, 16-20]: 



Vk(v) — A kz{ri) 



1-k 2 I" ^ 



z 2 {vi) 



At) z (if) 



■B k z(rj) 



"_df/_ 

z 2 (fj) 



z 2 (fj) 



(14) 



where Ak and Bk are fc-dependent constants that are de- 
termined by the Bunch-Davies initial condition (12) im- 
posed in the sub-Hubblc limit. As is well known, the first 
term involving Ak represents the growing mode, while the 
second containing B k corresponds to the decaying mode. 
However, it is important to realize that, at the time of 
preheating, the effective potential z" / z is no longer given 
by the slow roll expression (13). It is clear that the ef- 
fective potential will contain oscillatory functions and, 
hence, it can even possibly vanish. So, it is not a priori 
obvious that one can use the same approach as in the 
inflationary epoch and simply ignore the k 2 term in the 
differential equation (11) for arriving at the behavior of 
the super-Hubble modes. Moreover, it is known that, 
during the preheating phase, one has to deal with the 
resonant behavior exhibited by the equation of motion 
under certain conditions [29, 30, 77]. As a consequence, 
at this stage, it becomes necessary that we remain cau- 
tious and analyze equation (11) more carefully. 

In order to study the perturbations during the pre- 
heating phase, it proves to be more convenient to work 
in terms of cosmic time and use a new rescaled variable 
Vfc that is related to the Mukhanov-Sasaki variable as fol- 
lows: Vk = a 1 / 2 Vk- Then, one finds that Eq. (11) takes 
the form [29, 30] 
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dV 



EM 2 dej) 



V k = 0. (15) 



Recall that, in the quadratic potential of our interest, 
soon after inflation, the evolution of the scalar field <f>(t) 
is given by Eq. (4). Using this solution and its deriva- 
tive (5), it is then easy to show that, while the third, 
fourth and the fifth terms within the square brackets in 
the above differential equation decay as a -3 , the last term 
decays more slowly as it scales as a -3 / 2 . Upon retaining 
only the first, second and the last terms and neglecting 
the others, one arrives at an equation of the form 



d 2 V k 



da 2 



_±_ (^L) 3/2 cos (2cr + 2 A) 
mtp \ a J 



V k = 0, (16) 



where the new independent variable a is a dimensionless 
quantity which we have defined to be a = m t + 7r/4. We 
can rewrite the above equation as 



d 2 V k 
da 2 



[A k -2qcoB (2a + 2 A)] V fc = 0, (17) 



with Ak and q being given by 



<l 



= 1 



3/2 



(18) 
(19) 



where, as we mentioned, t e and a c denote the cosmic 
time and the scale factor when inflation ends. The above 
equation is similar in form to the Mathieu equation (see, 
for instance, Ref. [78]). The Mathieu equation pos- 
sesses unstable solutions that are known to grow rapidly 
when the values of the parameters fall in certain domains 
known as the resonant bands. As discussed in detail in 
Refs. [28, 30], since q <C 1 in the situation of our interest, 
one falls in the narrow resonance regime. In such a case, 
the first instability band is delineated by the condition 
1 — q < Ak < 1 + q, which turns out to be equivalent to 
the condition 



< - < V3Hm. 

a 



(20) 



It should be emphasized here that the time evolution of 
the quantities Ak and q are such that, once a mode has 
entered the resonance band, it remains inside it during 
the entire oscillatory phase. 

Note that, in Eq. (16), we can neglect the term involv- 
ing k 2 provided k 2 /(m 2 a 2 ) <C 1. This condition can be 
rewritten as 



k 



H 2 

— 5- < 1. 

m 



(21) 



On the other hand, the condition to fall in the first in- 
stability band, viz. Eq. (20), can be expressed as [30] 



k 



H 

3 m 



< 1. 



(22) 
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Given that, H < m immediately after inflation, it is ev- 
ident that the first of the above two conditions will be 
satisfied if the second is. In other words, being in the 
first instability band implies that one can indeed neglect 
the k 2 term in Eq. (16). But, clearly, this is completely 
equivalent to ignoring the k 2 term in the original equa- 
tion (f f ). Therefore, we can conclude that, provided we 
fall in the first instability band (which is the case for the 
range of modes and parameters of our interest), it is per- 
fectly valid to work with the super-Hubble solution (14) 
even during the preheating phase. 

The above conclusion can also supported by the fol- 
lowing arguments. As discussed in Ref. [30], in the first 
instability band, the Floquet index is given by /Lt = q/2. 
In such a case, the mode behaves as Vk cx e^ a . How- 
ever, in the situation of our interest, since we have a time 
dependent Floquet index, the corresponding solution can 
be written as 



Vfc cx exp 



/i der 



,3/2 



(23) 



which, in turn, implies that Vk = Vfc/a 1 / 2 cx a. Further, 
since, fk = v^/z and z cx a during preheating (i.e. if one 
makes use of the fact that ei = 3/2 on the average), we 
arrive at the result that fk remains a constant during this 
phase. This property is indeed the well known behavior 
one obtains if one simply retains the very first term of the 
growing mode in the super-Hubble solution (14). But, it 
should be realized that the above arguments also demon- 
strate another important point. Note that fk is a con- 
stant not only for modes on super-Hubble scales but, for 
all the modes (even those that remain in the sub-Hubblc 
domain), provided they fall in the resonance band. This 
behavior can possibly be attributed to the background. 
After all, it is common knowledge that the amplitude of 
the curvature perturbations remain constant on all scales 
in a matter dominated era [1-4, 6, 8-11, 16-20]. 

We have proven that, in the first instability band and 
on super-Hubble scales, the solution (14) is valid during 
preheating. Let us now analyze this solution in further 
detail. In particular, in order to check the extent of its 
validity, let us compare the analytical estimate for the 
curvature perturbation with the numerical solution. It is 



clear that the solution (14) leads to the following expres- 
sion for the growing mode: 



fh{v) - A * 



drj 

: 2 (ff) 



dfj z (fj) 



, (24) 



where we have retained the scale dependent correction 
for comparison with the numerical result. Since it is the 
contribution due to the growing mode that will prove 
to be dominant, we shall compare the behavior of the 
above solution during preheating with the corresponding 
numerical result. We shall carry out the comparison for 
a suitably small scale mode so that the second term in 
the above expression is not completely insignificant. 

We now need to evaluate the double integral in the 
above expression for fk during preheating. Let us write 



K(t) 
where 



" df > T(-\ 



z 2 (rj) 



dt 



J(t) 



dfj z 2 {fj) 



a(t)z 2 (t) 



dt 



J(t), (25) 



'tk z 2 (t)- (26) 
a(t) 



Upon making use of the matter dominated behavior of 
the scale factor and the expression (6) for e\, we find that 
the integral J(t) can be performed exactly. We obtain 
that 



J(t) = -M 2 pi a c t c 



6 ft 



5 U 



+ e 2iA {-2imt c )~°< a 7 ( -,-2imt 



5/3 



5/3 



+ e 



-2i A 



(2imt c ) 5/3 7 ( ^-,2imt ) +C 



, (27) 



where 7(6, x) is the incomplete Gamma function (see, for 
example, Refs. [79, 80]), while the quantity C is a di- 
mensionless constant of integration. Then, the relation 
between the incomplete Gamma and the Gamma func- 
tions allows us to express the function K(t) as follows: 



Kit) - 7^ 



dt 



4 a 2 J cos 2 (mt + A) 



6 ft 



1/3 



+ c~^ A (2imt c ) 



-5/3 



s If I 
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4/3 



-2i(mf+A) ( 2lmio )- 1 (|) +C 



4/3 



3tt /3 r* dit- 1 '* 



10 a 2 J cos 2 (mi + A) 



(28) 



In arriving at the final equality, we have used the asymp- totic property of the incomplete Gamma function [79, 80] 
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and have retained only the dominant term in inverse 
power of m t. The final expression above can be inte- 
grated by parts to arrive at 



K(t) 



3f 



4/3 



10 m al 



,.-1/3 



4/3 



10 m ai 



tan (m t + A) 
di i _4/3 tan (mi + A). (29) 



The second term containing the integral in this expres- 
sion is of the order of the other terms that we have al- 
ready neglected and, hence, it too can be ignored. As a 
result, the growing mode of the curvature perturbation 
can be written as 



At 



= A, 



3 kH 



2+4/3 



10 almt 1 ^ 



tan (m t + A) 



If k y H 
1- - — — tan (mt + A) 
5 V a tl J m 



(30) 



in perfect agreement with the result that has been ob- 
tained recently in the literature [32]. It is evident from 
the above expression that the evolution of the curvature 
perturbation will contain sharp spikes during preheat- 
ing, a feature that is clearly visible in Fig. 3 wherein we 
have plotted the above analytical expression as well as 
the corresponding numerical result (in this context, also 
see Fig. 4 in Ref. [30] where the spikes are also clearly 
visible). 

It is important that we make a couple of remarks con- 
cerning the appearance of the spikes in the evolution of 
the curvature perturbation. Firstly, as the spikes are en- 
countered both analytically and numerically, evidently, 
they are not artifacts of the adopted approach. Secondly, 
one may fear that the perturbation theory would break 
down as soon as one encounters a spike, which indicates 
a rather large value for the perturbation variable of inter- 
est. We believe that such issues could possibly be avoided 
when one couples the inflaton to radiation, as is needed 
to reheat the universe. 



vestigate the effects of preheating on the bi-spectrum. In 
this section, we shall quickly sketch the various contribu- 
tions to the bi-spectrum in the Maldacena formalism [37]. 

The scalar bi-spectrum B s (ki, k 2 , ^3) is defined in 
terms of the three point correlation function of the 
Fourier modes of the curvature perturbation 1Z as fol- 
lows [75, 76]: 



(K kl Kk 2 nk 3 ) = (2^) 3 B s (fei,fe 2 ,fe 3 ) 
x 6& (kx + k 2 + fe 3 ) • 



(31) 



10" 
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29.2 



29.4 29.6 
Number of e- 
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30.0 30.2 



FIG. 3: The behavior of the curvature perturbation during 
preheating. The blue curve denotes the numerical result, 
while the dashed red curve represents the analytical solu- 
tion (30). We have chosen to work with a very small scale 
mode k that leaves the Hubble radius at about two e-folds 
before the end of inflation. We have made use of the same 
value of A as in the previous two figures and we have fixed 
Ak [cf. Eq. (30)] by choosing it to be the numerical value 
of the curvature perturbation at a suitable time close to the 
end of inflation. It is clear that the agreement between the 
analytical and the numerical results is quite good. 



For convenience, we shall set 



III. THE SCALAR BI-SPECTRUM IN THE 
MALDACENA FORMALISM 

Now that we have understood the behavior of the large 
scale modes at the time of preheating, let us turn to in- 



G(kx,k 2 ,k 3 ) = (2Tr) 9 ^B s (k 1 ,k 2 ,k 3 ). (32) 

In the Maldacena formalism to calculate the bi- 
spectrum [37], the quantity G(ki, k 2 , k 3 ) can be ex- 
pressed as [38, 39, 43, 51, 52] 
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G(k 1} k 2 ,k 3 ) = ^2 G c (k 1 ,k 2} k 3 ) 



£7=1 



m 2 y 

PI / 4 



f kl (m) h 2 (m) h 3 (m) 5 C (fci , fe 2) fc 3 ) + far) ft 2 (m) ft 3 im) Q* (*i , to, fca) 



C=l 

■Gi(ki,k 2 ,kz) 



(33) 



where ?yf denotes the final time when the bi-spectrum is to be evaluated. The quantities Q c (fei, fc2, k 3 ) with C = (1, 6) 
are described by the integrals [38, 39, 43, 51, 52] 



Gi(kx,k 2 ,k 3 ) 
G2{ki,k 2 ,k 3 ) 

G 3 (k 1 ,k 2 ,k 3 ) 
G4{ki,k 2 ,k 3 ) 
Gs{ki,k- 2 ,k 3 ) 
G 6 {ki,k 2 ,k 3 ) 



2i dr) a 2 e\ (/^ /£ f' k \ + two permutations) , 

Jrii 

rm 

— 2i (fei • k 2 + two permutations) / dry a 2 t\ f kl f k2 f£ 3 , 

J rii 



rm 

2i di]a 2 el 
Jm 

rm 



kl u2 2 ) fki fk 3 fk 3 + fivc permutations 



rm 

/ d?y a 2 ex e' 2 (f k± f k2 f' k * 3 + two permutations) , 

J rii 

kx-k 2 



Arj a- e{ 



drj a 2 e{ 



k 2 



/fei fk 2 fk 3 + nve permutations 



k\ {k 2 ■ k 3 



h 2 h 2 



/fex /fe* 2 4*3 + tw ° permutations 



(34) 
(35) 

(36) 
(37) 
(38) 
(39) 



where rji denotes the time when the modes fk are well in- 
side the Hubble radius during inflation. The additional, 
seventh term Gr(ki,k 2 ,k 3 ) arises due to a field redefi- 
nition, and its contribution to G(ki,k 2 ,k 3 ) is found to 
be 



G 7 (fel,fe 2 ,fe3) 



£2(??f) 



two permutations 



(40) 



IV. THE CONTRIBUTIONS TO THE SCALAR 
BI-SPECTRUM DURING PREHEATING 

Our goal now is to determine the different contribu- 
tions to the bi-spectrum due to the epoch of preheat- 
ing for modes of cosmological interest. Since we know 
the behavior of the slow roll parameters and the large 
scale modes during preheating, it is clear from the above 
expressions that, it is simply a matter of being able to 
evaluate the integrals involved. As we shall illustrate, it 
turns out to be possible to actually evaluate these inte- 
grals explicitly and thereby arrive at the contributions to 
the bi-spectrum due to preheating. 

During preheating, for the large scale modes (i.e. for 
k <C aH), the contribution to fk is dominated by the 



constant growing mode [viz. the very first term in 
Eq. (14)] so that one has 



fk — A k 



(41) 



As the corresponding derivative trivially vanishes, at the 
same order in k, the leading non-zero contribution to the 
quantity f' k is determined by the decaying mode, which 
is given by 



fk — B k 



dr] 
z 2 (fj) 



and, hence, 



fk 



B, 



9 9 '■ 



(42) 



(43) 



where we have set B k = B k /(2 M 2 ). 

Let us now turn to discuss the various contributions to 
the bi-spectrum as the background and the modes evolve 
during preheating. 



A. The fourth and the seventh terms 

As we had mentioned in the introductory section, it 
has been realized that departures from slow roll during 
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inflation can lead to significant non-Gaussianities [44- 
52]. In such situations, it has been repeatedly noticed 
that it is the fourth term, viz. G4(fci, k 2 , k 3 ), that con- 
tributes the most to the bi-spectrum because it involves 
the derivative of the second slow roll parameter e 2 . Since 
e 2 grows extremely large during preheating [see Eq. (7) 
as well as Fig. 2], it is natural to expect that the fourth 
term will contribute significantly to the bi-spectrum at 
the time of preheating. So, we shall first investigate the 
contribution due to Gi(ki, k 2l k 3 ). As we shall see, it in- 
deed proves to be large during preheating. However, in- 
terestingly, as we shall illustrate, this large contribution 
is canceled by a similar contribution due to the seventh 
term G7(fci, k 2 , k 3 ) which also involves the second slow 
roll parameter e 2 . 

Upon using the supcr-Hubblc behavior (41) and (43) 
of the mode f k and its derivative in the expression (37), 
we obtain that 

£ 4 (fci, fc 2 , fc 3 ) ~ i (A* ki A* k2 Bl 3 + two permutations) 
rm 

x / dr,e' 2 , (44) 

J Tie 

where r/ c denotes the time at which inflation ends. The 
above expression can be trivially integrated to yield 

G<t{ki,k 2 ,k 3 ) — i (A% 1 A* k2 B k3 + two permutations) 
x [e 2 (r/ f )-e 2 (r/ e )], (45) 

so that the corresponding contribution to the bi- 
spectrum can be expressed as 

Gi(ki,k2,kg) - iM l M?7f)-e2(r?e)] 



\A kl \ 2 \A k2 \ 2 (A k3 B* k3 -A* k3 B k3 



two permutations 



(46) 



Since this expression is proportional to e 2 , it suggests 
that preheating may contribute substantially to the bi- 
spectrum. But, as we shall soon show, this large con- 
tribution is canceled by a similar contribution from the 
seventh term that arises due to the field redefinition 
[cf. Eq. (40)]. 
Now, consider the Wronskian 



w = hf:-fUL- 



(47) 



Upon using the equation of motion (9) for f k , one can 
show that, W = W/z 2 , where W is a constant. It is 
important to note that this result is valid on all scales, 
even in the sub-Hubble limit during inflation. In this 
limit, as we had mentioned, the modes v k satisfy the 
Bunch-Davies initial condition (12). On making use of 
this sub-Hubble behavior in the above definition of the 
Wronskian VV, one obtains that W = i. In the super- 
Hubble limit, we have, on using the corresponding solu- 
tion (41) and its derivative (43), 



W 



2 M 2 



Therefore, we obtain that 

(A k B* k -A* k B k ) 



2AP 



(49) 



and, hence, the expression (46) for G4(ki,k 2 ,k 3 ) simpli- 
fies to 



Gi(ki, k 2 ,k 3 



\ M?7f) - e 2 {r]e)} 
\A kl \ 2 \A k2 \ 2 



two permutations 



(50) 



Note that the first of these terms [involving e 2 (?yf)] exactly 
cancels the contribution Gi(k 1 , k 2 , k 3 ) [cf. Eq. (40)] that 
arises due to the field redefinition (with f k set to A k ), 
leaving behind only the contributions generated during 
inflation. 

Before we go on to discuss the behavior of the other 
contributions, we should emphasize here that the above 
result for the fourth and the seventh terms applies to all 
single field models. It is important to appreciate the fact 
that we have made no assumptions whatsoever about the 
inflationary potential in arriving at the above conclusion. 
However, one should keep in mind that, regarding its 
behavior near the minima, we have made use of the fact 
that the potential can be approximated by a parabola. 
Indeed, it is with this explicit form that we have been able 
to identify a solution to the Mukhanov-Sasaki equation 
that leads to a constant curvature perturbation. 



B. The second term 

We shall now turn to the second term G 2 (k 1 ,k 2 ,k 3 ), 
which in certain situations is known to be as large as the 
fourth term when there exist periods of fast roll [51, 52]. 
During preheating, we have 

G2(ki,k 2 ,k 3 ) = — 2i (fei ■ k 2 + two permutations) 

*A* kl A* k2 A* k3 I 2 (m,Vo), (51) 

where I 2 (rj{ 7 r] e ) denotes the integral 

f m 



rm 

I dria 2 el, 



(52) 



(A k B* k ~A* k B k ) 



(48) 



so that the corresponding contribution to the bi- 
spectrum is given by 

G 2 (ki, k 2 , k 3 ) = —2 i M pl (k± ■ k 2 + two permutations) 
x \A kl \ 2 \A k2 \ 2 \A k3 \ 2 
x [hiVfiVj-IZiVfiVe)] . (53) 

which identically vanishes since I 2 is real. Needless to 
add, this implies that the second term does not con- 
tribute to the bi-spectrum during preheating. Again we 



10 



should emphasize the fact that, as in the case of the 
fourth and the seventh terms, this result holds good for 
any inflationary model provided it can be approximated 
by a parabola in the vicinity of its minimum. 

In order to check that our assumptions and approxima- 
tions are indeed valid, let us now estimate the quantity 
02(feij &2, ^3) analytically during preheating for the case 



of the quadratic potential and compare with the corre- 
sponding numerical result. In such a case, the integral 
hitli, r] e ) can be carried out along similar lines to the in- 
tegral J(t) that we had evaluated earlier [cf. Eq. (27)]. 
We find that it can be expressed in terms of the incom- 
plete Gamma function 7(6, x) as follows: 



^(™g- 5/3 y (mt e 



,5/3 



e 5 (JV,-iV e )/2 _ l 



+ 4(-2i)- 5 / 3 c 24A 1 7 
+ 4(2i)- 5 / 3 e" 2lA 1 7 
+ (-4z)- 5 / 3 c 4iA (7 



-,-2 l mt c 3 W'-^)/ 2 
3 



7 ( -,-2imt 



3' 



2imt P e 3 W-^/ 2 



7 ( -,2imt ( 



7 ( -, -Aimtt 



(4i)- 5 / 3 e" 4lA <i7 



3 



-7 ( 2' 4im *< 



(54) 



On the other hand, had we ignored the oscillations during 
preheating, and assumed that the background behavior 
is exactly the same as in a matter dominated era, then, 
since, ((ei)) =3/2, the quantity ((hivf , Vc))) can be triv- 
ially evaluated to yield 



27 a e t e 
20 



3 5 (iV f -JVe)/2 



(55) 



We have plotted the quantity ^2(^1,^2,^3) in the equi- 
lateral limit, i.e. when k\ = k% = k^ — fc, corresponding 
to the analytical expressions (54) and (55) as well as the 
numerical result as a function of upper limit Nf during 
preheating in Fig. 4. The agreement between the analyt- 
ical and the numerical results is indeed striking. 




28.0 28.5 29.0 29.5 30.0 30.5 31.0 
Number of e— folds N 



C. The remaining terms 

In this sub-section, we shall analytically compute the 
contributions due to the remaining terms, viz. the first, 
third, fifth and the sixth. Notice that, the first term 
G\ (fci, /c2, /C3) and the third term Gs(ki, &2, ^3) involve 
the same integrals. Therefore, these two contributions 
to the bi-spectrum can be clubbed together. Similarly, 
the fifth and the sixth terms, viz. G^fci, &2, ^3) and 
Ge(ki,k2,k3), also contain integrals of the same type, 
and hence their contributions too can be combined. 

On making use of the super-Hubble behavior (41) 
and (43) of the mode and its derivative, we obtain 



FIG. 4: The behavior of the quantity ^2(^1, fe, ^3) in the 
equilateral limit, i.e. when k\ = — k$ — k, for a mode that 
leaves the Hubble radius at about 20 e-folds before the end 
of inflation. The blue curve represents the numerical result. 
The dashed red curve denotes the analytical result arrived 
at using the integral (54) and with the same choice of A as 
in the earlier figures. The dotted green curve corresponds 
to the integral (55) obtained when the oscillations have been 
ignored. As in the previous figure, the value of Ak has been 
fixed by choosing it to be the numerical value of the curvature 
perturbation on super-Hubble scales. Needless to add, the 
match between the analytical and the numerical results is 
excellent. 
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that, during preheating, 

Qi(ki,k 2l k 3 ) ~ 2 i (A* ki B^ 2 Bl 3 + two permutations) 
xha(m,Ve) (56) 



and 

G 3 {k 1 ,k 2 ,k 3 ) ~ -2i 



1 f e 2 2 J ^fei 5 fe2 ^3 



five permutations 



Il3(Vi,Vc)- (57) 



The quantity Ii 3 (j]f,rj e ) represents the integral 

t < \ r dr > 



(58) 



which can be trivially carried out during preheating to 
yield 



1 - e S(N f -N^)/2 



(59) 



From these results, we find that the contribution to the 
bi-spectrum due to the first and the third terms can be 
written as 



Gi(fci,fc 2 ,fc 3 ) + G 3 (fci,fc 2 ,fc 3 ) = 2iM? /i 3 (r?f,r/e) 



fci • fc 2 fei • fc 3 i 
1 — r^— I \A kl \ 



h 2 
h 2 



k 3 



{ A k 2 B* k2 A k;i B* k3 ~ A* k2 B k2 A* k;i B k3 ) + two permutations 



(60) 



Since the second term in the above expression for Ii 3 (rjf,r] e ) dies quickly with growing Nf, the corresponding contri- 
bution to the bi-spectrum proves to be negligible. 

The contributions due to the fifth and the sixth terms during preheating can be arrived at in a similar fashion. We 
obtain that 



G 5 (fci, fc 2 , fc 3 ) + G 6 (ki, fe 2 , fe 3 



iM 2 



fci • k 2 ki ■ k 3 k 2 (fe 2 • fe 3 ) 



k 2 



K 2 ""3 



I^J 2 (^ 2 5* 2 A k3 B* k3 



A* k2 B k2 A* k B k A + two permutations 



(61) 



with IsQ(r]f,r] e ) denoting the integral 



drj 



This integral too can be evaluated rather easily to arrive at the following expression: 



(62) 



3 ^ 



~HN t -N e )/2 cog 2 



heiVfiVe) = — cos (mt c + A)~e 

+ mt e cos(2A) |si(2mi e )-Si 2fflt c e 3(,v, ^=' /2 | 
+ mi c sin (2 A) jci(2mi e ) - Ci 2mt e e 3(Wf - JVe)/2 }), 



(63) 



where Si(a;) and Ci(x) are the sine and the cosine integral 
functions [79, 80]. And, had we ignored the oscillations, 



we would have arrived at 



3 1 
2al 



-3 (N t -N B )/2 



(64) 
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which is of the same order as /13 (r/f , r] e ), and hence com- 
pletely negligible as we had discussed. 



D. The contribution to / NL during preheating 

Let us now actually estimate the extent of the con- 
tribution to the non-Gaussianity parameter / NL during 
preheating. Since the contributions due to the combina- 
tion of the fourth plus the seventh and the second term 
completely vanish at late times, the non-zero contribu- 
tion to the bi-spectrum during preheating is determined 
by the first, third, fifth and the sixth terms. Note that, 
if one ignores the oscillations post-inflation, then one has 
keiViiVe) = 3 Iwivi 1 T /o)/2. In such a situation, we find 
that the non-trivial contributions lead to the following 
bi-spectrum in the simpler case of the equilateral limit: 



G oq (fc) 



69 i M 2 



I 13 ( m ,ri c )\A k \ 2 



(A 2 B* k 2 - Al 2 Bl) 



(65) 



In the equilateral limit, the non-Gaussianity parameter 
/nl is g iven by 



10 1 fc 6 G C q(fc) 

"9"(2^r7 P s 2 ( fc ) ' 



(66) 



where V s {k) is the power spectrum defined in Eq. (8). 
Upon making use of the fact that f k ~ Ay, at late times, 
we then obtain the contribution to during preheating 
to be 



115 i M 2 



48 

AlBl 2 



k 2 m 



\A k 



(67) 



In order to explicitly calculate the parameter / NL , we 
need to first specify the inflationary scenario. We shall 
choose to work with power law inflation because it per- 
mits exact calculations, and it can also mimic slow roll 
inflation. During power law inflation, the scale factor can 
be written as a(rj) = a\ (t?/?7i) where a\ and 771 are 
constants, while /? is a free index. It is useful to note that, 



in such a case, the first slow roll parameter is a constant 
and is given by ei = {fi + 2)/(/3 + 1). The current con- 
straints on the scalar spectral index suggest that /3 < — 2, 
which implies that the corresponding scale factor is close 
to that of de Sitter. 

In power law inflation, the exact solution to Eq. (11) 
can be expressed in terms of the Bessel function J v (x) as 
follows: 



v k (r)) = V Z &0 [C k M-kr)) + D k J_ u (-kri)}, (68) 

where v = ((3 + 1/2), and the quantities C k and D k 
are constants that are determined by the initial condi- 
tions. Upon demanding that the above solution satisfies 
the Bunch-Davics initial condition (12), one obtains that 



k 2 cos (tt /?) ' 



(69) 
(70) 



One can confirm that the super-Hubble limit of the so- 
lution (68) above indeed reproduces Eq. (14) exactly. 
Moreover, the limit also allows us to arrive at the con- 
stants Ak and Bk, which are found to be 



B k = - 



Ck, 



r(0 + 3/2) V277aiM pl 
(2/3 + l)2' 3 + 1 /2 ^2T iai M p 



r(-/? + i/2) 

x (-k m )- D k . 



(71) 



(72) 



Then, upon inserting the above expressions for the quan- 
tities Ak and Bk in Eq. (67), we find that 

/cq (fc ) = Hhl T ^((3 + l) 2^+1 (2/3 + l) 2 
Jn ^ ' 288tt V 2/ 



x sin(2 7r^) |/3 + l 



-2(/3+i) 



-3 (Ni-N B )/2 



-(2/3+1) 



(73) 



This expression can also be rewritten in terms of the pa- 
rameters describing the post-inflationary evolution. We 
obtain that 



-3 (JV f -JV e )/2 



TT 2 9* 

30 



-1/4 



1/41 -(2/^+1) 
1/4 Peri 



a H 



-(2/3+1) 



(74) 



where denotes the effective number of relativistic de- grees of freedom at reheating, Ti-h the reheating tern- 
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perature and z eq the redshift at the epoch of equality. 
Also, p cr i, ao and H represent the critical energy den- 
sity, the scale factor and the Hubble parameter today, 
respectively. The above expression is mainly determined 
by the ratio p~J^ /T r h- For a model with /3 ~ — 2 and 
a reheating temperature of T r h — 10 10 GeV, one obtains 
that / NL ~ 10 -60 for the modes of cosmological interest 
(i.e. for k such that k/ao ~ Hq), a value which is com- 
pletely unobservable. This confirms and quantifies our 
result that, in the case of single field inflation, the epoch 
of preheating does not alter the amplitude of the scalar 
bi-spectrum generated during inflation [81]. However, it 
is worthwhile to add that, while the amplitude of the 
above non-Gaussianity parameter / NL is small, it seems 
to be strongly scale dependent. 

We believe that a couple of points require further em- 
phasis at this stage of our discussion. Recall that, to fall 
within the first instability band during preheating, the 
modes need to satisfy the condition (20). But, in order 
to neglect the term involving k 2 in the differential equa- 
tion (16), the modes of interest are actually required to 
satisfy the condition (21). As we have emphasized ear- 
lier, evidently, the condition (21) will be easily satisfied 
by the large scale modes that already lie within the in- 
stability band and thereby satisfying the condition (22). 
Therefore, it is important to appreciate the fact that the 
conclusions we have arrived at above apply to all cosmo- 
logically relevant scales. 

V. DISCUSSION 

In this work, we have analyzed the effects of preheat- 
ing on the primordial bi-spectrum in inflationary mod- 
els involving a single canonical scalar field. We have il- 
lustrated that, certain contributions to the bi-spectrum, 
such as those due to the combination of the fourth and 
the seventh terms and that due to the second term, van- 
ish identically at late times. Further, assuming the infla- 
tionary potential to be quadratic around its minimum, we 
have shown that the remaining contributions to the bi- 
spectrum are completely insignificant during the epoch 
of preheating when the scalar field is oscillating at the 
bottom of the potential immediately after inflation. It is 
important to appreciate the fact that the results we have 
arrived at apply to any single field inflationary potential 
that has a parabolic shape near the minimum. 

A couple of other points also need to be stressed re- 
garding the conclusions we have arrived at. The results 
we have obtained supplement the earlier results wherein 
it has been shown that the power spectrum generated 
during inflation remains unaffected during the epoch 
of preheating (see, Ref. [29]; in this context, also see 
Ref. [82]). Moreover, our results are in support of earlier 



discussions which had pointed to the fact that the contri- 
butions to the correlation functions at late times will be 
insignificant if the interaction terms in the actions at the 
cubic and the higher orders depend on either a time or a 
spatial derivative of the curvature perturbation [83-87]. 

Broadly, our effort needs to be extended in two dif- 
ferent directions. Firstly, it is important to confirm that 
the conclusions we have arrived at hold true for potentials 
which behave differently, say, quartically, near the min- 
ima. Further, the exercise needs to be repeated for mod- 
els involving non-canonical scalar fields [88-90]. In this 
context, it is worth mentioning that the generalization of 
the conserved quantity TZk in the Dirac-Born-Infeld case 
has been shown to stay constant in amplitude on scales 
larger than the sonic horizon, a property which allows us 
to propagate the spectrum from horizon exit till the be- 
ginning of the radiation dominated era [89] . Secondly, as 
we had mentioned, preheating is followed by an epoch of 
reheating when the energy from the inflaton is expected 
to be transferred to radiation. It will be interesting to 
examine the evolution of the bi-spectrum during reheat- 
ing. However, in order to achieve reheating, the scalar 
field needs to be coupled to radiation. It is clear that 
the formalism for evaluating the bi-spectrum involving 
just the inflaton is required to be extended to a situation 
wherein radiation too is present and is also coupled to 
the scalar field. 

However, possibly, the most interesting direction 
opened up by our work concerns multi-field inflation and 
associated non-Gaussianities [91-96]. Unlike single field 
models wherein the curvature perturbation associated 
with the large scale modes is conserved at late times, such 
a behavior is not necessarily true in multi-field inflation. 
When many fields are present, the entropy (i.e. the iso- 
curvature) fluctuations can cause the evolution of curva- 
ture perturbations even on super-Hubble scales. Further, 
in the case of multi-field inflation, it is known that the 
two-point correlation function can be affected by met- 
ric preheating [77]. In other words, the power spectrum 
calculated at the end of multi-field inflation is not nec- 
essarily the power spectrum observed in, say, the CMB 
data because the post-inflationary dynamics (that is to 
say preheating instabilities) can modify it. In the same 
manner, it is tempting to conjecture that the scalar bi- 
spectrum calculated at the end of multi-field inflation will 
not necessarily be the same as the one observed in the 
data [97]. In particular, the so-called consistency rela- 
tions [43] , which relate the three-point and the four-point 
correlation functions (or, equivalently, the corresponding 
dimensionless non-Gaussianity parameters / NL and t nl ) 
might receive corrections in a multi-field context due to 
metric preheating. We are currently investigating such 
issues. 
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